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Adiabatic state engineering is a powerful technique in quantum information and quantum control. However, 
its performance is limited by the adiabatic theorem of quantum mechanics. In this scenario, shortcuts to adi- 
abaticity, such as provided by the superadiabatic theory, constitute a valuable tool to speed up the adiabatic 
quantum behavior. Here, we propose a superadiabatic route to implement universal quantum computation. Our 
method is based on the realization of piecewise controlled superadiabatic evolutions. Remarkably, they can be 
obtained by simple time-independent counter-diabatic Hamiltonians. In particular, we discuss the implementa¬ 
tion of fast rotation gates and arbitrary n-qubit controlled gates, which can be used to design different sets of 
universal quantum gates. Concerning the energy cost of the superadiabatic implementation, we show that it is 
dictated by the quantum speed limit, providing an upper bound for the corresponding adiabatic counterparts. 


I. INTRODUCTION 

Quantum adiabatic processes are a powerful strategy to im¬ 
plement quantum state engineering, which aims at manipu¬ 
lating a quantum system to attain a target state at a designed 
time T. In the adiabatic scenario, the quantum system evolves 
under a sufficiently slowly-varying Hamiltonian, which pre¬ 
vents changes in the populations of the energy eigenlevels. In 
particular, if the system is prepared in an eigenstate |n(0)) of 
the Hamiltonian H at a time f = 0, it will evolve to the cor¬ 
responding instantaneous eigenstate \n{t)) at later times. This 
transitionless evolution is ensured by the adiabatic theorem, 
which is one of the oldest and most explored tools in quan¬ 
tum mechanics UHll . The huge amount of applications of 
the adiabatic behavior has motivated renewed interest in the 
adiabatic theorem, which has implied in its rigorous formu¬ 
lation H HTOll as well as in new bounds for adiabaticity IfTTI - 
E] . In quantum information processing, the adiabatic the¬ 
orem is the basis for the methodology of adiabatic quantum 
computation (AQC) 11141 , which has been originally proposed 
as an approach for the solution of hard combinatorial search 
problems. More generally, AQC has been proved to be uni¬ 
versal for quantum computing, being equivalent to the stan¬ 
dard circuit model of quantum computation up to polynomial 
resource-overhead na . Moreover, it is a physically appeal¬ 
ing approach, with a number of experimental implementations 
in distinct architectures, e.g., nuclear magnetic resonance M- 
[18] , ion traps IIT9ll , and superconducting flux quantum bits 
(qubits) through the D-Wave quantum annealer Il20l4^ . 

Recently, the circuit model has been directly connected 
with AQC via hybrid approaches 112^ l24ll . Then, an adia¬ 
batic circuit can be designed based on the adiabatic realiza¬ 
tion of quantum gates, which allows for the translation of the 
quantum circuit to the AQC framework with no further re¬ 
sources required. In particular, it is possible to implement 
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universal sets of quantum gates through controlled adiabatic 
evolutions (CAE) Il24ll . In turn, CAE are used to perform one- 
qubit and two-qubit gates, allowing for universality through 
the set of one-qubit rotations joint with an entangling two- 
qubit gate Il25ll26l . However, since these processes are ruled 
by the adiabatic approximation, it turns out that each gate of 
the adiabatic circuit will be implemented within some fixed 
probability (for a finite evolution time). Moreover, the time for 
performing each individual gate will be bounded from below 
by the adiabatic time condition Bno] . Eor a recent analysis 
on adiabatic control of quantum gates and its corresponding 
non-adiabatic errors, see Ref. ||271 . 

In order to resolve the limitations of adiabaticity in the 
hybrid model, we propose here a general shortcut to CAE 
through simple time-independent counter-diabatic assistant 
Hamiltonians within the framework of the superadiabatic the¬ 
ory GHHM] . The physical resources spent by this strategy 
will be governed by the quantum circuit complexity, but no 
adiabatic constraint will be required in the individual imple¬ 
mentation of the quantum gates. Moreover, the gates will be 
deterministically implemented with probability one as long 
as decoherence effects can be avoided. In particular, we dis¬ 
cuss the realization of rotation gates and arbitrary n-qubit con¬ 
trolled gates, which can be used to design different sets of 
universal quantum gates. This analog approach allows for fast 
implementation of individual gates, whose time consumption 
is only dictated by the quantum speed limit (QSL) (for closed 
systems, see Refs. iMia). Indeed, the time demanded for 
each gate will imply in an energy cost, which increases with 
the speed of the evolution. In this context, by analyzing the 
energy-time complementarity, we will show that the QSL pro¬ 
vides an energy cost for superadiabatic evolutions that upper 
bounds the cost of adiabatic implementations. 

II. ADIABATIC QUANTUM CIRCUITS 

Let us begin by discussing the design of adiabatic quantum 
circuits as introduced by Hen ll24i through the implementation 
of quantum gates via CAE. 
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A. Controlled adiabatic evolution 

In order to define quantum gates through CAE, we will 
introduce a discrete bipartite system associated with a 

Hilbert space 'Hs ® 'Hji- The system is composed by a 

target subsystem S and an auxiliary subsystem whose in¬ 
dividual Hilbert spaces and 'Kji have dimensions ds and 
dj{, respectively. The dynamics of will be governed by a 
Hamiltonian in the form Il24ll 

it) [l ® + g (f) Pk ® ], (1) 

where /(O) = g(l) = 1, g(0) = /(I) = 0, and {Pk} denotes 
a complete set of orthogonal projectors over S, so that they 
satisfy PkPm - ^kmPk YikPk - 1 ■ Alternatively, we can 
write Eq. ([T]| as 

= (2) 

with Ht (t) - g (t) +f(t) Pi'-'’'' denoting a Hamiltonian that 

acts on Suppose now that we prepare the system SJl in the 
initial state = \^)®\eb), where |i/r) is an arbitrary state of 
S and \sb) is the (non-degenerate) ground state of Then 
I'i'inir) is the ground state of the initial Hamiltonian By 

applying the adiabatic theorem Ena, a sufficiently slowing- 
varying evolution of H(t) will drive the system (up to a phase) 
to the final state 

Y,,Pkm»\sk}, (3) 

where lej.) is the ground state of m ■ 

1. Single-qubit unitaries and controlled two-qubit gates 

We can perform a single-qubit unitary transformation 
through a general rotation of an angle tp around a direction 
h on the Bloch sphere. In this direction, we begin by prepar¬ 
ing the system SlPl, taken here as two qubits, in the initial 
state \^) ® |0), where {|0), |1)) are the computational states of 
the auxiliary system Then, we let the system adiabatically 
evolve driven by the Hamiltonian Il24l 

H{s) - |n+)<n+| ®Ho(s) -H |n_) <n_| ®(s) , (4) 

where Hq(s) and H^{s) are adiabatically-evolved Hamiltoni¬ 

ans, whose effect will be restricted to the respective sub¬ 
spaces of the projectors |n±){n±| = ^ (!+«■(?), where 

« is a unitary vector on the Bloch sphere associated with 
S and (T = (o-x,(Ty,cr^^, with {cr,) denoting the set of 
Pauli matrices. The Hamiltonians are taken as H^(s) - 
-tdh jcTj cos(0O'S) + sinjffo'S) cos ^ -t- cr^ sin j, where ^ 6 
{0, (p), cob sets the energy scale (co > 0), 0o is a constant pa¬ 
rameter, and i = f/r denotes a dimensionless (parametrized) 
time, with t the total time of evolution. Note then that 

Ho (s) = -cob {cTj cos (0os) + CTji sin (0qs)] , 

H^ (s) = -tob {cTj cos(0os) + sin(0os) [cr^^ cos cp -i- cry sin (^j|. 

( 5 ) 


By writing the initial state of SJl as = |iA>®|0> = 

(a |n+) -H j3 |n_))® |0), where |i/f) is an arbitrary (not necessarily 
known) qubit state, the final state follows from Eq. ([^, 

i.e. = a\n+} ® +/3\n^} ® |e“). Note that is 

the ground state of H^(s), reading = cos (0O'S/2) |0) -i- 
e'^ sin ( 6*0 v/2) 11), with ^ 6 {0, tp}. An equivalent form of writ¬ 
ing \'¥final) is given by 

I'y final) = cos (6»o/2) (a\h+) -Hy0|«-»® |0) 

+ sin (0o/2) [a |n+) + \n.}) ® |1). (6) 

Hence, we have a probabilistic implementation of the rotated 
target state a\n+) ■+■ e“'’l3\hP) for an arbitrary angle cp around 
an arbitrarily chosen axis n, with probability sin^ (0o/2). In 
particular, this probability approximates to one by taking 0o ~ 

7T. 

In order to perform controlled rotations of a qubit by an 
angle cp around a direction n, the starting point will be to take 
the subsystem as a two-qubit system and keeping ^ as a 
single auxiliary qubit. The Hamiltonian is now chosen to be 

H(s) - (|0,n+><0,n+| H- |0,n_)(0,n_| H- |1, n+) (1, n+|) ® Hq (i) 
+ |l,«-)<l,n-|®//^(v), (7) 

which will govern the evolution of the initial composite state 
IT'mif) = |iA)®|0>^ with \ip) = Q;|0,n+)-Hy6|0,n_)-Hy|l,n+) -I- 
d|l,n_). From Eq. the final state of the subsystem S 
in the limit 0o ^ ^ is now the controlled rotated vector 
_ Q,|Q^ |1,«_). By com¬ 

bining controlled rotations with the single-qubit unitaries dis¬ 
cussed above, it is possible to design universal sets of quantum 
gates through an adiabatic implementation. 

III. RESULTS 

In this Section we present the main results of this work. We 
start by generalizing the adiabatic implementation of quan¬ 
tum gates proposed in Ref. M for n-qubit controlled gates. 
Even though n-qubit controlled gates can be decomposable 
into one and two-qubit gates (see, e.g. Refs. ||2511^ ). this 
implementation implies in an extended class of adiabatic uni¬ 
versal gates, e.g. the set {Toffoli, Hadamard) llJTlI^ . Then, 
we will derive the main result of this work, which is a short¬ 
cut for general adiabatic circuits through constant counter- 
diabatic Hamiltonians, which implies in the possibility of fast 
analog implementations of quantum circuits. Moreover, we 
will present an analysis of the quantum speed limit in the con¬ 
text of the energetic cost of the superadiabatic circuit. 

A. Adiabatic «-controlled gates 

In order to implement n-controlled gates, we define the sub¬ 
system S as an («-i- l)-qubit system, with the first n qubits used 
as the control register and the last qubit taken as the target reg¬ 
ister. For the auxiliary system we keep it as a single qubit. 
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Then, we take the initial state as I'T;,,,-,) = \ij/) ® |0), with the 
subsystem S described by 

l<A) = Jki,...,k„,e\k\,---,kn,h^), (8) 

where yki....,k„,e are complex amplitudes, ki e {0,1), e - {+), 
and h is an arbitrary axis in the Bloch sphere. Here we have 
written the target qubit in the basis {|n±)), leaving the remain¬ 
ing qubits of S in the computational basis. For simplicity, we 
will write the states in its decimal representation, i.e. 


B. Shortcut to adiabaticity via counter-diabatic driving 

Let us introduce now a shortcut to general CAE through the 
superadiabatic approach. This will allow for fast piecewise 
implementation of quantum gates, whose evolution time will 
not be constrained by the adiabatic theorem. We begin by 
defining the evolution operator 





ni 


A^-I 

EE ym,e \m, Hf) , 

m=Q €=± 


where N = 2". Then, we let the system evolve driven by the 
Hamiltonian 


H(s)^ 


N-2 

^ ^ |n, «£>(«, I + |A^- l,n+)<A^- 1,/ 

n=0 €=± 

+ \N -l,n_}{N - l,n_|®//^(v). 


® Ho (s) 

(9) 


which leads an initial state |'T(0)) = |n(0)) into an evolved 
state I'T (0) given by 

where |n) = \n (f)) are the eigenvectors of the adiabatic Hamil¬ 
tonian. Note that this evolution mimics the adiabatic behavior. 
The Hamiltonian that guides the evolution of the system is the 
superadiabatic Hamiltonian, which reads 

HsA(t) = H(t) + HcD(t), (15) 


We note that the rotation of the target qubit is expected to be 
applied if the state of the control system is \N - 1). Then, if 
the Hamiltonian is sufficiently slowly-varying so that we can 
apply the adiabatic theorem, the system will achieve the final 
state 

N-2 

final) = \m,K) ®\sq) + yN-\,+ W - l,n+)® Ifig) 

m=0 e=± 

-i-yw_i_|A-l,n_)® |e^), (10) 

where |£^)is defined as = cos (0o/2) |0)-i-e'^ sin(0o/2) |1) 
(^ 6 {0,0)). An equivalent form of writing Eq. is 

\^fina) = COS (0o/2) ^ ® |0) + sin (0o/2) ® |1> , (11) 

with 

£=± 

+ |A - 1) ® (rw-i,+ l«r) + yN-f-e'^ \n-}) . (12) 

Thus, by performing a measurement over the auxiliary qubit, 
we find the rotated state with probability sin^ iOo/2). 

As in the previous case of a rotation controlled by one qubit, 
this probability can be enhanced to one in the limit 6o —> tt. 
Indeed, this state implies in a rotation of the target qubit in 
S conditioned by the state |A-1) s |l---l)of the control 
register in S. An application of this scheme is the adiabatic 
implementation of the Toffoli gate, which constitutes an uni¬ 
tary operation implementing an X gate over the target qubit 
if all control qubits are in the state 1, with no effect if any 
qubit of the control register is in the state 0. This can be easily 
achieved here by a rotation of an angle n around of the di¬ 
rection X, therefore choosing (p - ti and |n±) = |+), with |±) 
denoting the eigenstates of cr^. 


where the additional term Hcd (0 is the counter-diabatic 
Hamiltonian 

Hcd it) = ib ^ i\d,n) {n\ H- {d,n\n) l«> <«l) ■ (16) 

n 

Therefore, a superadiabatic implementation of a dynamical 
evolution involves the knowledge of the eigenstates of the adi¬ 
abatic Hamiltonian H (f), which limits the direct application of 
the superadiabatic approach in quantum computation. Eor in¬ 
stance, by adopting the original AQC approach lfT4ll , supera¬ 
diabatic implementations seem prohibitive, since the whole 
set of eigenlevels of a many-body Hamiltonian is required. 
In a similar context, counter-diabatic driving protocols based 
on realizable settings have been investigated for assisted evo¬ 
lutions in quantum critical phenomena Il39l - l4n . Here, as 
we shall see, the superadiabatic implementation of univer¬ 
sal quantum circuits in the hybrid approach can be promptly 
achieved, since we deal with the eigenspectrum of piecewise 
Hamiltonians, which act over a few qubits. It is then appealing 
to formulate a superadiabatic theory to CAE and then to spec¬ 
ify it to the implementation of universal sets of quantum gates. 
Let us begin by establishing the complete set of eigenstates of 
the Hamiltonian H (t) provided by Eq. Q. To this end, con¬ 
sider the eigenvalue equation to each Hamiltonian Hk (f) given 
by 

Hk(t)\4(t))^Ei(t)\s[(t)), (17) 

with / = 1, • • • , djf. By defining the projectors Pk in Eq. Q as 
Pk - |d/t) {Ak\, with {Ak\Ak') - dkk' and k - I,-■ ■ ,ds ,'nq can 
write the complete set of eigenstates of H (t) as 

|r;(o) = i^z)®H(f)), (18) 

such that H(t)\y‘i(t)) = Ej (0 |y) (O)- Indeed, from Eq. (|^, 
we have H(t)\y\(t)) = [E/t® (0] [l-I/) ® |e;(f))] = 
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E'l (t) |yj (f)). Note that each projector is associated with 
a Hamiltonian Hk. For instance, for the adiabatic imple¬ 
mentation of n-controlled gates, we have defined the Hamil¬ 
tonian H in Eq. 0 by linking the set - 

0, - ■ ■ ,N - 2), \N - l,n+}{N - l,n+|) with Hq and by link¬ 
ing the remaining projector \N - l,nP){N - \,nP\ with H^. 
The next step is to obtain the counter-diabatic Hamiltonian 
Hcd (0 that implements the shortcut to the adiabatic evolu¬ 
tion of //(f). In this direction, we use the eigenstates of //(f) 
as given by Eq. ([T^. Then, we get 





with s jyl (f)^. Therefore 


quantum dynamics as provided by the quantum speed limit 
(QSL). We will consider a closed quantum system evolving 
between arbitrary pure states |'T (0)) and |'T (t)). Since the 
evolution is driven by a time-dependent superadiabatic Hamil¬ 
tonian HsAit), we will take the generalized Margolus-Levitin 
bound ll^ derived by Deffner and Lutz llT5l , which reads 

^ JcosX(T'(0),T'(t))-1| 
r>n - - -, (24) 

tjT 

where X ('T (0), T* (t)) = arccos [K'Pol'I' (t))|] is the Bures 
metric for pure states and 

Er^- f f/fK'T(0)|//sA(OI'P(f)>|. (25) 

T Jo 


Hcd (0 = ifi ^ 

i 





= ;^jp,®//™(f)], 



where s |ej (f)^ and is the counter-diabatic Hamilto¬ 
nian to be associated with the piecewise adiabatic contribution 
Hi (f) acting over subsystem J/ 1 , which reads 


Eor superadiabatic evolutions, the initial state IT* (0)) = l7o(0)) 
evolves to |'T (f)) = lyoCO)^ where lyoCO) denotes the instanta¬ 
neous ground state of the adiabatic Hamiltonian //(f). By us¬ 
ing the parametrized time s - f/r, we can show from Eqs. ( [24| ) 
and ( |25| ) that the total time r that mimics the adiabatic evo¬ 
lution within the superadiabatic approach can be reduced to 
an arbitrary small value. More specifically, the addition of a 
counter-diabatic Hamiltonian implies into the QSL bound 


//f" (f) = m ^. [|5,£') (e;| -H (d,si 



( 21 ) 


Hence, from Eq. ( fTS] !, we can implement the shortcut dynam¬ 
ics through the superadiabatic Hamiltonian 


HsA{t)^Y^Jk®Hl^{t), ( 22 ) 


where it) = Hk (f) -H (f) is the piecewise supera¬ 
diabatic Hamiltonian. Note that the cost of performing su¬ 
peradiabatic evolutions requires the knowledge of the eigen¬ 
values and eigenstates of Hi (f). Eor the implementation of 
general «-controlled gates, this is a Hamiltonian acting over 
a single qubit, which is independent of the circuit complex¬ 
ity. Moreover, we can show that, for an arbitrary n-controlled 
quantum gate, the counter-diabatic Hamiltonians //^^ (f € 
{ 0 , 0 )) associated with shortcuts to adiabatic evolutions driven 
by H^ (s) = —u)ii {cTj cos( 0 os) + sin( 0 os) [cTr cos ^ + cTy sin^jj, 
with s - t/T, are time-independent operators given by 

//^® - \^y cos cTx sin . (23) 

Eq. ( [23] l shows that the implementation of the shortcut can 
be achieved with a very simple assistant Hamiltonian in the 
quantum dynamics. Its proof is provided in Section Methods. 


C. Quantum speed limit 

It is expected that the shortcut via a counter-diabatic Hamil¬ 
tonian is faster than the evolution via an adiabatic Hamilto¬ 
nian, but how much faster can it be? To answer this ques¬ 
tion, we shall take a lower bound to the time evolution in 


f/Tw-HA- > ^|cosX('1'(0),'1'(t))-1|, (26) 

with 77 > 0 and ^ ^ |cos £. ('T ( 0 ), T* (t)) - 11, as shown in 

Section Methods. Therefore, the QSL bound reduces to 

TU>0 (V l'l'(O)) , |'T(t)>), (27) 


with T and w defined by the superadiabatic Hamiltonian 
HsaJ)- This means that the superadiabatic implementation 
is compatible with an arbitrary reduction of the total time t, 
which holds independently of the boundary states IT* ( 0 )) and 
IT* (t)). Naturally, a higher energetic cost is expected to be in¬ 
volved for a smaller evolution time t. In particular, saturation 
of Eq. ( [27] i is achieved for either t —» 0 or oj ^ 0, with both 
cases implying in to —» 0. Note that this limit is forbidden 
in the adiabatic regime for finite co, since the energy gap is 
proportional to ho, which implies in an adiabatic time of the 
order Tad 1/w", with n e ll^ 14^171 14211 . Hence, Eq. (271 

leads to a flexible running time in a superadiabatic implemen¬ 
tation, only limited by the energy-time complementarity. 


D. The Energetic Cost 

Let us show now that time and energy are complementary 
resources in superadiabatic implementations of quantum evo¬ 
lutions. We shall define the energetic cost associated with a 
superadiabatic Hamiltonian through 

E(t)=- r WHsAimdt, ( 28 ) 

T Jo 

with Hsa (0 given by Eq. (j2^ and the norm provided by the 
Hilbert-Schmidt norm ||A|| = sjTr [A'l'A]. Since Hsa (0 is 
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Hermitian, we can write 


= ^^r[HHt) + Hl^{t)\dt. 


(29) 


To derive Eq. (291, we have used that Tr(|//(f), //c£>(0}) = 0. 
This can be obtained by computing the trace in the eigenbasis 
of Hit) and noticing that the expectation value of Hcoit) taken 
in an eigenstate of //(f) vanishes, i.e. {y'iit)\ Hcoit) \y\{t)') — 0 . 
In particular, let us define the energetic cost to the adiabatic 
Hamiltonian as 


^Tr[H2it)]dt. 

Then, it follows that the energetic cost 2 (t) in superadiabatic 
evolutions supersedes the energetic cost Eq (t) for a corre¬ 
sponding adiabatic physical process. In order to evaluate E (t) 
we adopt the basis of eigenstates of the adiabatic Hamiltonian 
Hit). By using Eq. ( fTS] !, this yields 

^ W = ? £ a/IIm + *V," (Opr, (30) 

where E'J' it) are the energies of the adiabatic Hamiltonian 
Hi it) and 


eT (O = {9,7? (0 \dt7T (O) - |(rr (0 \d,yl' (0)f ■ (31) 


In order to analyze the energetic cost as provided by Eq. 
for superadiabatic qubit rotation gates, we set ds - dj{ - 2 
and E'” is) = (-l)™^' hto (V/). Moreover, by using Eq. (18 1 , 


we 


obtain i^fit) = {d,£fit)\d,s^^it)) - \{s'^it)\d,e1'it))\ , 
which leads to jj."' is) - 0g/4T^ [See Eqs. (34i and ( [35] l in 
Section Mef/io(/s]. Hence 


E(t) = 2. 1 + 




4(tw)2 




(32) 


We illustrate the behavior of E (t) in Eig. where it is ap¬ 
parent that the energetic cost increases inversely proportional 
to the total time of evolution. In particular, note also that, for 
a fixed energetic cost, the optimal choice Oo tt requires a 
longer evolution. This is because of the fact that, in this case, 
the final state associated with the auxiliary qubit is orthogonal 
to its initial state, so it is farther in the Bloch sphere. In the 
more general case of controlled gates, the analysis is similar 
as in the case of single-qubit gates. However, we must take 
into account the number of projectors composing the set {Pk}- 


More specifically, the sum over / in Eq. (301 shall run over 
1 to 4, which is the number of projectors over the subsystem 
S. Thus we can show that energetic cost to implement 
controlled gates is E*''^ = V 2 E. 



FIG. 1. Energetic cost in unities of few as a function of wt for differ¬ 
ent values of 8 q . 


IV. METHODS 

A. Time-independent connter-diabatic Hamiltonians for 
u-controlled gates 

Let us explicitly design here the superadiabatic implemen¬ 
tation of controlled evolutions for piecewise Hamiltonians 
H^is) as provided by Eqs. 0 . To this end, consider the eigen¬ 
value equation 


H^is)^is)) = E‘^is)\e^is)), (33) 


with ^ = (0, (fi}, where 

14 (i)) = - sin ^ |0) -i-e'^cos ^ |1>, (34) 

|4 (■*)) = ^ I0) + sin ^ 11>. (35) 


Erom Eq. a, it follows that the eigenstates for the adi¬ 
abatic Hamiltonian His) governing the composite system 
SJl are given by the sets (s)^ = |n+) ® |eg (s)^| and 


(i)^ = |«-) ® |sjj (■s))| associated with the 


set of eigen¬ 


values j/ig (i)| and \e'^ (i)|, respectively. By evaluating the 
eigenvalues of //q is) and H^ is), we obtain that their spectra 
are equal, being provided by E^ — E"^ - ±a)h. Thus, His) 


exhibits doubly degenerate levels, with ||yg (s)^ , (■*))! 

{Ito (^))>| 7<p ('*^)| associated with levels E^ - wfe and E = 
-wife, respectively. By using now Eqs. (34i and (35 i, we ob¬ 
tain (dts'i I e'S — 0, for any i — {+) and ^ = (0,0). Then, 
from Eq. (20 1 , we obtain that the counter-diabatic hamilto- 
nian is H^^ is) - ih Zi/=|±) (s)^ (e'^ (s)|, which leads to 

the time-independent counter-diabatic Hamiltonian given by 
Eq. ([2^. The extension to the case of n-controlled gates can 
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be achieved as follows. From Eq. (18 i, the eigenstates of H (s) 
read 


that |^fj<ro( 0 )lro(s)>l > l<i .9 Kro ( 0 ) Iro (s)>ll (see proof in 


Iri 


0(iV-l) * 



Ref 1351), which yields 


)) = N,Mf)® 14 (■*))’ 



)) = IN - 1,77+)® 14 (i)), 

xis) > mis) >h f ds 

ds Kro (0) Iro ( 5 ))! 

)) = IN - 1,77_)® |e^(5)). 

Jo 

pi 



where m - (0, • • • , - 2), e,k - {+) 

By computing the eigenvalues of 
tain that the spectrum of H{s) is 

with 


and ^ 
H(s), 


= ( 0 , 0 ). 

we ob- 


[\rol ’ |K^v-i) (^)) > |r^^v-i) 


(2A^)-degenerate, 


and 


{|4„ (^)) ’ |r^ov-i) (^)) ’ |r^(v-i) 


| 7 ^(^ associated with the 

levels - ojfi and - -a>H, respectively. By using these 
results into Eq. (20i, we obtain that the counter-diabatic 


piecewise Hamiltonian (s) is given by Eq. (|23|l. Hence, 


the implementation any n-controlled gate is achieved through 
a time-independent counter-diabatic Hamiltonian. 


> fi 


f 


ds id, Kro (0) Iro (s)>l) 
= ^IKro(0)lro(i))l-i|, 


(38) 


where we have used the inequality dt\f(t)\ > dtf(t)\. 
Erom the dehnition of the Bures metric, we have 
Kro(0)lro(l))l = cosi:(yo(0),ro(l))- Hence, ^(i) > 


h |cos X. (i/' (0), >f/ (t)) - 11, which implies into Eq. (27 1 . 


V. DISCUSSION AND CONCLUSION 

B. Quantum speed limit for super adiabatic evolutions 


Let us apply here the QSL bound to superadiabatic evolu¬ 
tions. By using the fact than the I*? it)} evolves in the ground 
state lyo (0) of H (f) and that Hgyi (t) is given by Eq. ( [T5] l, we 
have 


Er 



dt\{yoiO)\HsA (0 Iro (0)1, 

dt |£o (0 <ro (0) Iro (0) + <ro (O) \Hcd (0 Iro (0)1, 


where Edit) is the instantaneous ground state energy of El{t). 
Now we use Eq. (16 1 and the inequality dt \f (t) + g (x)| < 

£ dt \f (01 + £ dt |g (01, which yields 


Er<- f dtlEo(t)(ro(0)lro(t)}l + - f dt\{yomd,jom 
T Jo T Jo 

-H- r c/f Kro (OIJrro(0><ro(0) Iro (0)1 ■ (36) 

T Jo 

By using the parametrized time s = tfr, we obtain 


Er < 771 (s) -H 


(mis) + 7 / 3 ( 5 )) 

T 


where the parameters rjiis) (/ = 1,2,3) are given by rf\{s) - 
£ ds \Eq is) (yo (0) Iro (s))l, m(s) = ds Kro (0) IJ^ro (s))! 

and 7 / 3 (s) = ^ Kro (s) I J^ro (■«)) <ro (0) Iro (s))l- Since 

the ground state energy for the adiabatic Hamiltonian //(s) 
in the case of 77 -controlled gates is Eo(s) = —a>h [see 
Eqs. 0 and 0 ]^ we write 7/i(s) = corjis), with 7/(s) = 

t/s Kro (0) Iro (s))!- Moreover, we dehne;^^(s) s 7 / 2 ( 5 ) -H 
7 / 3 ( 5 ). Then 


rj(s) ojT + xis) > /7|cosX(ro(0),ro(l)) - 1| ■ (37) 


Let us now analise the term 7 ( 5 ). Eirst, note that 
(ro (0) IJsro (s)) = «ro (0) Iro (5))). Then, we use 


We have proposed a scheme for implementing universal 
sets of quantum gates within the superadiabatic approach. In 
particular, we have shown that this can be achieved by ap¬ 
plying a time-independent counter-diabatic Hamiltonian in the 
auxiliary qubit to induce fast controlled evolutions. Remark¬ 
ably, this Hamiltonian is universal, holding both for perform¬ 
ing single-qubit and 77-controlled qubit gates. Therefore, a 
shortcut to the adiabatic implementation of quantum gates can 
be achieved through a rather simple mechanism. In particu¬ 
lar, different sets of universal quantum gates can be designed 
by using essentially the same counter-diabatic Hamiltonian. 
Moreover, we have shown that the flexibility of the evolution 
time in a superadiabatic dynamics can be directly traced back 
from the QSL bound. In this context, the running time is only 
constrained by the energetic cost of the superadiabatic im¬ 
plementation, within a time-energy complementarity relation¬ 
ship. Implications of the superadiabatic approach under de¬ 
coherence and a fault-tolerance analysis of superadiabatic cir¬ 
cuits are further challenges of immediate interest. In a quan¬ 
tum open-systems scenario, there is a compromise between 
the time required by adiabaticity and the decoherence time 
of the quantum device. Therefore, the superadiabatic imple¬ 
mentation may provide a direction to obtain an optimal run¬ 
ning time for the quantum algorithm while keeping an inher¬ 
ent protection against decoherence. In turn, a basis for such 
development may be provided by the generalization of the su¬ 
peradiabatic theory for the context of open systems P3l - l46l . 
Concerning error-protection, it may also be fruitful the com¬ 
parison of our approach with non-adiabatic holonomic quan¬ 
tum computation, where non-adiabatic geometric phases are 
used to perform universal quantum gates (see, e.g. recent pro¬ 
posals in Refs. Il47ll48l ). Moreover, the behavior of correla¬ 
tions such as entanglement may also be an additional relevant 
resource for superadiabaticity applied to quantum computa¬ 
tion. These investigations as well as experimental proposals 
for superadiabatic circuits are left for future research. 
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